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Abstract 

An accelerator experiment was performed using a low-energy antiproton beam 
to measure antiproton detection efficiency of BESS, a balloon-borne spectrome- 
ter with a superconducting solenoid. Measured efficiencies showed good agreement 
with calculated ones derived from the BESS Monte Carlo simulation based on 
GEANt/gheisha. With detailed verification of the BESS simulation, the relative 
systematic error of detection efficiency derived from the BESS simulation has been 
determined to be ±5 %, compared with the previous estimation of ±15 % which 
was the dominant uncertainty for measurements of cosmic-ray antiproton flux. 
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1 Introduction 



The BESS spectrometer, shown in Fig. 1, was designed [1,2] and developed [3,4,5,6] 
as a high-resolution balloon-borne spectrometer with the capability to search 
for rare cosmic-rays and provide various precision measurements of cosmic- 
ray primaries. In spite of many advantages in cylindrical configuration, the 
adoption of solenoidal magnets could not be realized in previous balloon-borne 
spectrometers because of much unavoidable material in particle passage. How- 
ever, a thin superconducting solenoid developed at KEK [2,4] enabled us to 
uniquely adopt this concentric configuration. A uniform magnetic field of 1 
Tesla is produced by a thin superconducting coil [4], through which particles 
can pass without too many interactions. The magnetic-field region is filled 
with a tracking detectors (JET/IDCs), resulting in an acceptance of 0.3 m 2 sr. 
Tracking is performed by fitting up to 28 hit points in the drift chambers. The 
upper and lower scintillator hodoscopes (TOF) [6] provide time-of-flight and 
two dE/dx measurements. The instrument also incorporates a threshold-type 
Cherenkov counter [5] with a silica-aerogel radiator to distinguish high-energy 
p's from electron and muon backgrounds. 
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Fig. 1. Cross-sectional view of the BESS detector showing a p event. 

Since 1993 through 2000, seven balloon flights have been successfully carried 
out and more than 2 x 10 3 antiprotons (p's) have been unambiguously detected 
[7,8,9,10,11]. This has allowed measuring the energy spectrum of cosmic p's 
and investigating their origin. To investigate the origin of low energy p's more 
sensitively, it is inevitably important to reduce systematic errors as well as 
statistical errors of the resultant spectrum. The dominant source of systematic 
error in the low energy region (< 1 GeV) is the uncertainty in p interaction 
losses in the instrument. In the previous analyses, we surmised ±15 % relative 
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error to the detection efficiency (e) , which is defined as 

£ -^obs/-^inc; 

where iV inc is the number of incidence within the acceptance of the detector, 
and iV bs is the number of identified particles. Since A^" inc cannot be derived 
directly from flight data, the detection efficiency is evaluated using the Monte 
Carlo simulation (BESS MC) [12] based on the geant/gheisha code [13,14]. 
The BESS MC incorporates detailed material and detector descriptions such 
that realistic detector performance is obtained. The original gheisha code 
was modified so that experimental data of p-nuclei cross sections are repro- 
duced [15,16]. However, it is difficult to estimate the systematic error due to 
interaction losses because of uncertainties in secondary multiplicity, angular 
distribution, and detector response; accordingly, detection efficiency must be 
directly measured and must precisely verified to reduce systematic error in the 
low-energy region p flux. 

Considering this, we performed an accelerator beam experiment at the KEE- 
PS K2 beam line using a low-energy p and proton (p) beam. The objectives 
of the beam experiment were as follows: 

(1) directly measure detection efficiencies for p's and p's; 

(2) examine the BESS MC simulation; 

(3) reduce systematic error in detection efficiency especially for p's. 

Although the BESS detector has been successively upgraded since the first 
successful flight in 1993, the basic detector concept is the same, i.e., (i) large- 
acceptance cylindrical configuration with solenoidal magnet, and (ii) mass 
identification using the rigidity and velocity measurements. Accordingly, the 
presented beam results can be applied to both past and future detectors. 

Section 2 summarizes experimental setup, after which Section 3 describes inci- 
dent beam identification which determines N inc . Measurement of the detection 
efficiency is then presented in Section 4, with a study of simulation results be- 
ing discussed in Section 5. Finally, Section 6 summarizes the main results and 
provides conclusions. 



2 Experimental Setup 

The BESS beam experiment was performed in February 1999 at the KEE- 
PS K2 beam line which is equipped with an electro-static separator [17] to 
enrich low energy p's. Figure 2 shows a schematic view of the experimental 
setup at the down stream of the K2 beam line. D2 is a dipole magnet, Q6 and 
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Q7 are quadrapole magnets. They are used for beam transport and focusing. 
KURAMA is another dipole magnet to analyze momentum of the incident 
particles. The BESS detector was rotated circumferentially 70° in r-0 plane 
(a plane perpendicular to the axis of the solenoid) such that it was suitably 
positioned in the beam line for proper (from top to bottom) beam incidence. 
To identify incident particles and to reject interacted events, we placed four 
trigger counters (T1-T4), two drift chambers (DC1 and DC2), and an aerogel 
Cherenkov counter (AC) in the beam line (for details, see figure caption). 
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Fig. 2. Experimental setup. Tl - T4 are 10 mm thick plastic scintillators, having 
time resolutions of 30 - 40 ps. DCs use a mixture of 50 - 50 % Ar-ethane gas with 
a full drift length of 1 cm. Each DC has three planes of horizontal wires and three 
planes of vertical wires. The spatial resolution is 150 /im/wire. AC uses 8-cm thick 
aerogel radiator with a refractive index of 1.03. 



Data from the BESS detector and beam line detectors are collected using 
the BESS data acquisition system. To obtain N inc , instead of using the BESS 
standard trigger (coincidence between upper and lower TOF counters), the 
trigger was generated by Tl & T2 & T3 & T4 & AC. Note that AC remarkably 
improves the p/ir ratio. 

Data were collected for three different detector configurations to represent the 
typical incidence of cosmic-ray particles in terms of the amount of material 
and penetrated region. Three configurations, CFG1-3, are shown in Fig. 3 
together with typical p trajectories. The kinetic energy of incident particles at 
the BESS top of instrument (£"toi) ranges from 0.1 to 1 GeV (0.4 to 1 GeV) 
for p's Qo's). The BESS detector cannot be rotated more than 70 ° due to a 
constraint of the internal structure of the liquid helium storage. Therefore, low 
energy p's were out of the BESS acceptance region in this beam experiment 
due to the opposite deflection. 
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Fig. 3. Cross-sectional views of the BESS detector with incident p beams for three 
configurations (CFG 1-3). Using CFG1 as a reference, CFG2 is the case in which the 
incident beam has an angle of cos# = 0.915, where 6 is defined as the angle between 
the beam and y-axis in the y-z plane. CFG3 is the case in which the incident beam 
passes through a different region of the central tracking system in r-<f> plane. The 
distance of incident position between CFG1 and CFG3 is 153 mm along the x-axis. 

3 Analysis 

3.1 Beam identification 

To determine iVi nc the beam line detectors must unambiguously identify the 
incident particle and precisely determine incident energy, position, and angle. 
The incident kinetic energy, £7toi> was derived from /3t4-T3, where (3, the 
velocity of incident particle, was obtained from a pair of timing measurements. 
Note that (i) less accurate energy resolution was obtained from the deflection 
measurement due to the multiple scattering and limited track length, and (ii) 
it is not appropriate to derive the incident energy from the /3t4-ti, etc. due to 
the energy losses in the beam line in spite of the much better 1/(3 resolution. 
The accuracy of energy determination was 1 % in the very low energy region 
around 0.2 GeV. Around 1 GeV, energy determination is obtained with 4 % 
accuracy due to the constant 1//3t4-T3 resolution. Systematic error of the 
absolute energy is estimated to be ± 1 % due to calibration of the time of 
flight measurements and energy losses in the beam line. The beam trajectory 
obtained by DC1 and DC2 is extrapolated to the BESS detector taking into 
account the fringing field of the solenoid. The accuracy of the incident position 
and angle around 1 GeV were 2.5 mm and 2 mrad (rms deviation), respectively. 

To identify p's and p's, we required that the energy losses in Tl - T4 and 
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the velocities measured by various combination of trigger counters are con- 
sistent with p's. Figure 4 shows examples of these cuts for p's of -Etoi ~ 
1 GeV. Since 1/(5 distribution shows a clear separation between p, kaon, and 
pion/muon/electron particles which allows incident beam particles to be un- 
ambiguously identified. To ensure that incident particles arrive at the top of 
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Fig. 4. 1/(3 and dE/dx cut by trigger counters for p's of .Etoi ~ 1-0 GeV. 1//9t3-tij 
1//?T4-T3, and d£'/dxT4 are shown in (a), (b), and (c), respectively. The hatched 
histogram indicates events survived after applying the p beam selection. 

the BESS instrument, we further required (i) at least one hit in the upper 
TOF hodoscope and (ii) beam trajectory agreement with the upper TOF hit 
position in both r-0 and y-z planes, where the hit point in TOF is obtained 
from the segmentation of hodoscope in <fi direction and timing difference of 
PMT at the both end in z direction. These two requirements guaranteed that 
the incident particles pass through T4 and AC without large angle scattering 
or interaction, resulting in the precise determination of Ai nc . 

In order to compare the beam data and BESS MC results, a MC data set was 
generated as follows. To estimate interaction and energy losses in T4 and AC, 
which were located just upstream of BESS, they were described in the BESS 
MC. Input kinematics of beam particles was obtained from beam data event 
by event. This allows comparisons between beam data and MC results under 
the same conditions. Data sets are referred to as BEAM and MC samples, 
respectively. 



3.2 Detection efficiency 



In order to select non-interacting p's, we applied the same cuts as those used 
in the standard p analysis for the flight data [10]: 

(1) select events with a single downward-traveling particle fully contained in 
the fiducial region of the tracking detectors; 
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(2) require only one (one or two) hit in the upper (lower) TOF hodoscopes 

[18]; 

(3) require the TOF hit position consistent with the extrapolation of the 
trajectory determined by JET and IDCs (not by beam line trackers); 

(4) require that dE/dx measurements of upper and lower TOF are consistent 
with those of p's; 

(5) require a quality track in terms of fitting chi-squares, number of hits used 
in fitting, and etc.; 

(6) require that the mass derived from momentum and velocity measure- 
ments is consistent with p mass. 

Cuts (l)-(4) reject most of the events with interactions. Distributions of the 
cut parameters are provided in Ref. [19]. Cut (5) assures the quality of momen- 
tum and velocity measurements. Cut (6) is to identify p's. Selection criteria 
for p's are identical to the above except for the sign of momentum. The ef- 
ficiencies of cuts (l)-(4) are different for p's and p's because of the different 
interaction cross sections. However, the efficiencies of cuts (5) and (6) are the 
same for p's and p's due to the symmetrical configuration of the detector. 
Therefore, we express the detection efficiency (e) as the product of the non- 
interacting efficiency (£ n0 n-int) associated to the cuts (l)-(4) and the quality 
and particle-identification efficiency (£q_id) associated to cuts (5) and (6): 

£ = £non-int " £Q-ID- 

Since £ n0 n-int relies on the MC calculation and £q-id can be estimated by using 
the p sample in the flight data, we should verify the former efficiency by the 
beam experiment. 

3.3 Beam-related corrections and systematic errors 

The following beam related corrections are applied to £ non -int and the system- 
atic errors are estimated. 

(1) Beam dump effect: Annihilation of p's at the beam dump located 3 m 
downstream of the BESS detector produces secondary particles, some of 
which generate delayed hits in the TOF counters. After eliminating the 
beam dump effect by removing the delayed hits from the hits in the upper 
and lower TOF hodoscopes, the remaining correction was negligible. The 
systematic uncertainty was estimated to be 0.002. 

(2) Accidental tracks: Multi-track events without a vertex in the detector 
were considered to be a beam-related accidental track and the track is 
removed from the event. However, an accidental track would be mis- 
identified as an interaction if the accidental particle is close to the incident 
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p. The remaining correction for mis-identification was O.OOf - 0.005 and 
the systematic uncertainty of the £ n on-mt was estimated to be 0.001. The 
correction and its error for p's were negligible because of the smaller 
incident beam intensity. 

(3) Interaction at upper TOF: The use of upper TOF information to identify 
the non-interacting beam entering the BESS detector rejects a small frac- 
tion of events with interaction inside BESS, which underestimates N[ nc 
and thus overestimates e non _ int . The correction for this effect was esti- 
mated by varying the cut parameters of the selections and was estimated 
to be —0.01 0.005. The rms deviation of the difference in £ n0 n-int be- 
tween BEAM and MC results was taken as the systematic uncertainty 
which was evaluated to be 0.005. 

(4) Beam energy: Since the systematic uncertainty of the beam energy is 
±1 %, the corresponding uncertainty in £ n0 n-int is negligible for ener- 
gies above 0.16 GeV. However, it is difficult to reliably estimate the 
corresponding systematic errors in £ non -int derived from BEAM below 
0.16 GeV, where £ n0 n-int changes rapidly with energy. 

In total, beam related systematic uncertainty was estimated to be less than 
0.01 for energies up to 1 GeV except for very low energy region where the 
efficiency rapidly changes with energy. 



4 Results 

Figure 5 shows £ non -int for p's and p's for each configuration derived from 
BEAM and MC samples. Reflecting the energy dependence of cross sections, 
£ n0 n-int f° r p' s gradually decreases with decreasing energy, while e non -int for 
p's shows little dependence on energy. Among BEAM data, the differences 
between CFG1-3 in e n0 n-int are 7, 5, and 3 % at energies of 0.2, 0.4, and 
1 GeV, respectively. Although large discrepancies exist below 0.16 GeV where 
^non-int rapidly drops due to stopping of incident p's in the instrument, they 
can be explained by the systematic uncertainty of 1 % in the beam energy 
determination. 

The relative differences in £ non -int between BEAM and MC samples ((Ae/s) non -[ 
(emc — ^bcamV^bcam) are shown in Fig. 6 from 0.16 to 1.0 GeV for p's and 0.4 
to 1.0 GeV for p's, where beam related systematic errors were kept low and 
well estimated. As shown, (Ae/e) non „ int was kept within ±5 % forp's, demon- 
strating good agreement between BEAM and MC samples. Above 0.5 GeV 
small but significant discrepancies (2 %) in (A£/£) non _i nt are present, while 
statistical accuracy of BEAM data limited the discussion in the low energy 
region below 0.3 GeV. For p's, (Ae/e) non _i nt was kept within ±2 %. Moreover, 
good reproducibility of (A£/£) non _ int are obtained between different configu- 
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Fig. 5. Direct measurement of £ non -int fork's andp's in (a) CFGI, (b) CFG2, and (c) 
CFG3. The error bars include both statistical and systematic uncertainties, while 
they are too small to see in this figure (see Fig. 6, instead). 
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Fig. 6. Relative difference in e n0 n-int of MC from BEAM for (a) p's and (b) p's. 

rations, i.e., the maximum deviation of Ae non _ int between all configurations 
are close to the statistical accuracy of each data point for all energies (mostly 
~ 0.005). Note that CFG1-3 were selected to represent the typical incidence 
of cosmic-ray particles, and the properties of non-interacting particles show 
very good agreement between BEAM and MC samples. 

As a result, these results clearly demonstrates that above 0.16 GeV (for p's 
above 0.4 GeV) systematic error in £ non -int is determined to be within ±5 % 
(±2 %). Since systematic error in £q_id was well estimated to be within ±1 % 
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using unbiased p sample in flight data, systematic error in £ non -int can be 
considered as that in detection efficiency. 



5 Discussion 



In a detailed study of the interaction processes in the BESS MC [19], we found 
no significant difference in £ non -int f° r different hadronic packages, gheisha 
and fluka [20] (implemented in GEANT3, denoted as gfluka) if we assume 
the same total inelastic cross section. This means that the multiplicity and 
angular distribution of secondary particles are less important as far as £ n on-mt 
is concerned. Any interactions in the BESS detector are well discriminated 
by a combination of various and independent measurements. The only factors 
relevant to the efficiency are essentially the total inelastic cross section and the 
material amount. We studied dE/dx of stopping cosmic rays using the flight 
data and compared it with the BESS MC simulation. As a result we confirmed 
that the amount of material was correctly implemented in the BESS MC with 
an accuracy of ±0.5 g/cm 2 . 

Since ~l/3 of the incident particles interact in the instrument, the ±5 % 
agreement between the measured and calculated efficiencies means that the 
previous measurements of p-nuclei inelastic cross sections were verified at least 
within an accuracy of about ±15 %. 

Figure 7 shows the inelastic and elastic cross sections for p to aluminum used 
in the BESS MC (gheisha(bess), solid curve), as well as those of the orig- 
inal GHEISHA code, FLUKA code, and experimental data [15]. The data of 
Kuzichev et al. only represent annihilation cross section. Assuming that the 
non-annihilation inelastic cross section is about 10 % of the total cross section 
as assumed in the gheisha code, we modified the data of Kuzichev et al. and 
refit all the data below 15 GeV, indicated by the dash-dotted curve in Fig. 7 
(gheisha (new)). With this new cross section, £ non -mt was reproduced with 
better accuracy up to 1 GeV. As a result, the systematic error in detection 
efficiency can be further reduced to ±2 % from 0.3 to 1 GeV for p's. 

Since we do not have the measured data of £ non -int above 1 GeV, the systematic 
error in this energy region was estimated from the relative difference in £ n on-mt 
among several hadronic packages (gheisha(new), gheisha(bess), gfluka). 
The systematic error in detection efficiency above 1 GeV was estimated to be 
±5 %. 



10 



.2 3000 - 



p + Aluminum 



Inelastic (upper) 
Elastic (lower) 



GHEISHA 
FLUKA 
DATA FIT 
DATA FIT (NEW) 




100 



■ Nakamura et al 
» Kuzichev et al 
' Gorin et al 
o Allaby et al 
n Afonasiev et al 
4 Carrol et al 



Elastic 



10 10 

Kinetic Energy (GeV) 



Fig. 7. Inelastic and elastic cross sections for p to aluminum. Histograms are from the 
GHEISHA and fluka codes and data points [15] from Nakamura et al. (filled squares); 
Kuzichev et al. (filled triangles); Gorin et al. (inverted triangles); Allaby et al. (open 
circles); Afonas'ev et al. (open squares); and Carroll et al. (open triangles). The solid 
curves are fit with these data points and used in the BESS MC. The dash-dotted 
curve is a new fit that better reproduces the results of the beam experiment. 



6 Conclusion 



Through an accelerator beam experiment, measurements of the detection ef- 
ficiency for p's and p's in BESS were performed below 1 GeV. The measure- 
ments and following studies indicated that the detection efficiency was mainly 
determined from material amount in the instrument and cross sections with 
nuclei, essentially regardless of the simulation details such as multiplicity and 
angular distribution of secondaries. As a result, we confirmed that the relative 
systematic error in detection efficiency derived from the simulation was kept 
below 5% and 2% for p's and p's, respectively. It can be applied to system- 
atic errors in cosmic-ray p measurements [8,9,10,11], as well as future BESS 
experiments including high statistics long-duration flights[21] providing the 
instrumental features of BESS are maintained. Moreover, we found the to- 
tal advantage of cylindrical configuration with superconducting solenoid by 
eliminating possible disadvantage in determining detection efficiencies accom- 
panied with 30-40 % p interaction losses in the instrument. By increasing the 
reliability of the cosmic-ray p spectrum, these results will enable us to carry 
out the most sensitive-ever investigation on the origin of cosmic-ray p's. 
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